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ABSTRACT 

We investigate signatures of population III (PopIII) stars in the metal-enriched environment 
of GRBs originating from population II-I (PopII/I) stars by using abundance ratios derived 
from numerical simulations that follow stellar evolution and chemical enrichment. We find 
that at z > 10 more than 10% of PopII/I GRBs explode in a medium previously enriched by 
PopIII stars (we refer to them as GRBII—HII). Although the formation of GRBII—HII is more 
frequent than that of pristine PopIII GRBs (GRBIIIs), we find that the expected GRBII—HII 
observed rate is comparable to that of GRBIIIs, due to the usually larger luminosities of these 
latter. GRBII—>-111 events take place preferentially in small proto-galaxies with stellar masses 
M* ~ 10 4 - 5 — 1O 7 M 0 , star formation rates SFR ~ 10 -3 — lO _1 M 0 /yr and metallicities 
Z ~ 10 -4 — 10 -2 Zq. On the other hand, galaxies with Z < 10 -28 Z 0 are dominated by 
metal enrichment from PopIII stars and should preferentially host GRBII—>-111. Hence, mea¬ 
sured GRB metal content below this limit could represent a strong evidence of enrichment by 
pristine stellar populations. We discuss how to discriminate PopIII metal enrichment on the 
basis of various abundance ratios observable in the spectra of GRBs’ afterglows. By employ¬ 
ing such anal ysis, we conclude t hat the currently kno wn candidates at redshift z ~ 6 - i.e. 
GRB 050904 ( iKawai et al.lf2006l) and GRB 130606A dCastro-Tirado et al.ll2013l) - are likely 
not originated in environments pre-enriched by PopIII sta rs. Abundance measurements_for 
GRBs at 2 ~ 5 - such as GRB 100219A dThone et alJl2013l) and GRB 111008A dSnarre et alJ 
12014 - are still poor to draw definitive conclusions, although their hosts seem to be dominated 
by PopII/I pollution and do not show evident signatures of massive PopIII pre-enrichment. 

Key words: cosmology: theory - early Universe - gamma-ray burst: general - galaxies: high 
redshift 


1 INTRODUCTION 

Primordial epochs host the formation of first stars and galaxies, 
and the consequent production of heavy elements in the Uni¬ 
verse. The basic understanding of the topic relies on cosmic 
molecule formation and gas cooling in growing dark-matter po¬ 
tential wells. Features of primordial structures are still debated, 
tho ugh, and their observational signatures under inve stigation (see 
e.g. lCiardi & Ferrai32005i ; lBromm & Yoshidalfeoi ll for recent re¬ 
views). Early star formation episodes are expected to occur at red¬ 
shift {z) larger than ~ 10, when the Universe was only half a 
Gyr old. In these ages, star formation is led by pristine gas col¬ 
lapse that drives the birth of the so-called population III (PopIII) 
stars. PopIII stars evolve and die ejecting the heavy elements they 
have synthesized in their cores. These events pollute the surround¬ 
ing medium and change its chemical composition. Thus, the fol¬ 
lowing generations of population II-I (PopII/I) stars will be born 
in a medium that has been pre-enriched of metals. The transition 
from the primordial PopIII to the subsequent PopII/I regime is im¬ 


portant to understand early cosmic phases, although the detailed 


(e.g. Bromm & Loeb 2003); Schneider et al. 

2003|; Tornatore et al. 

2007b; Maio et al. 20id; Peon et al] 2014 

; Yajima & Khochfat 

2014|), For instance, PopIII stars have been predicted to have both 


very large and small masses ( e.g. lSchneider et all2002l : lciark et alJ 
1 201 11 : IStacv & Brommll20f4l) . while PopII/I stars have theoretical 
metal yields that often are affected by stellar evolution modeling. 
The cosmological transition from one regime to the other is ex¬ 
pected to take place at very early times, with a PopIII contribu¬ 
tion to the cosmic st ar formation rate ( SFR) dropping d ramatically 
below 2 ~ 10 (e.g. iMaio et al]|20ld : IWise et aul2012h . The very 
first star formation sites are supposed to have dark-matter masses 
as small as ~ 2 x 10 6 M@, to be extremely bursty and to feature 
specific SFRs up t o ~ 10 — 100 Gyr" 1 i Salyaterraetal J [20 1 3l : 
iBiffi & Maid 1201 31: Ide Souza et aljlloll IWise et al.ll20l4 . Such 
objects are enriched to average metallicities of ~ 10 -2 Zq by 
a > 9. 
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Direct observations of primordial times are quite problematic, be¬ 
cause early (proto-)galaxies are usually faint and difficult to detect. 
Chemical imprints of PopIII stars may be retrieved from metal- 
poor stars of the Galactic halo dKarlsson et all2008l : lRollinde et alj 
1 20091) and low-mass halos at high redshift h osting damped 
Ly-al pha absorbers (D LAs)dT renti et al] 120091 : iKulkami et aid 
120131 : iMaio et all 12013l : IKulkami et al.l 2014h. but their features 


are s t ill largely deba t ed 1 Tumlinsonetaf 20041; DaijmeetalJ 


20041: lErni et al.l 120061: iTumlinsonl 

Trend & Shullll20ld:lKomiva et al.f 


200?riKobavashi et al.l 1201 ll: 

201 Oil. 


A powerful tool to study this epoch is represented by gamma-ray 
bursts (GRBs), true cosmic li ghthouses that have b e en observed 
up to 2^8 — 9 (s ee e.g . Isalvaterra et aid 120091 : iTanvir et all 
120091 : ICucchiara et aid 1201 ID and can be employed for several 
purposes. In particular, high-redshift GRBs can provide funda¬ 
mental information about the early stages of structure formation 
and the properties of their own hosti ng galaxies, such as (see 
McQuinn et al] |2009l : [Amati et a l~]|2013l for a m ore complete list) : 
metallicity and dust content 1 Sava2lioetalj2005|; Nuza_et al. 200 1 
Icampisi et ID l201ll : iMannucci et al.l 1201 ll: lElliott et aid boll): 
neutr a l-hydrogen fraction 1 Gallerani^t_a^]2(ffl8|^ ]Nagamine et all 


120081: iMcOuinn et alj|2008l: iRobertson & Ellisll2012h : local 


tnter- 


galacdc radiat i on fie ld _dlnoue_et alj_j20 1 d) : stellar populations 


(Camnisi et al. 

2011 

Salvaterra et al. 

2013 


Tomartjd] 1201 ll: Ide Souza et al. 1 201 ll ; 


early cosmic magnetic fields dlakahasM_etal][201jJ); dark-matter 


models dde Souza et al.ll2013l:jMaio & Viell 20141) : and primordial 


non-Gaussianities (Maio et aU l2012l) .' 


In the following, we will focus on the metal content of first 
GRB host sites in order to explore their gas chemical ratios and 
to constrain the corresponding stellar populations. This is inter¬ 
esting for observers that try to infer the original stellar popula¬ 
tion from the species detected in the spectra of GRB afterglows. 
This study will be performed by relying on published N-body hy- 
drodynamical chemistry simulations including atomic and molec¬ 
ular non-equilibrium cooling, star formation, feedback mecha¬ 
nisms, stellar evolution, metal spreading according to yields and 
lifetimes of different populations (PopIII and PopII/I) for sev- 
eral chemical_species jMaio et aDbOOlEoToi: ICampisi et al.l201 ll: 

ISalvaterra et al] 120131) . As the time evolution of the individual 

heavy elements is consistently followed within the simulation run, 
it will be possible to properly compute expected chemical ratios for 
PopIII and PopII/I GRB progenitors and to assess PopII/I GRBs 
exploding in environments possibly pre-enriched by former PopIII 
events. 

Throughout this work a standard ACDM cosmological model with 
total matter density parameter Do, m = 0.3, dark energy den¬ 
sity parameter Do ,a = 0.7, baryonic matter density parameter 
Do ,6 = 0.04, Hubble constant in units of 100 km s _1 Mpc -1 
h = 0.7, spectral normalization as = 0.9 and primordial spec¬ 
tral index n = 1 is assumed. 

The paper is organized as follows: in Section[2]we describe the sim¬ 
ulations employed in this work; in Section[3]we discuss the results 
in term of rate of GRBs, metal signatures and properties of the host 
galaxies; in Section[4]we give our Conclusions. 


2 SIMULATIONS AND METAL ABUNDANCES 

In this paper we use the simulations presented in IMaio et aP 
] 2010 l) . of which we describe the essential features, while we re¬ 
fer the reader to the original paper for more details. The same 


X 

logio(Afx/A H )g 

X 

logio(Afx/A H )g 

C 

-3.57 

Si 

-4.49 

N 

-4.17 

S 

-4.86 

O 

-3.31 

Ca 

-5.66 

Mg 

-4.4 

Fe 

-4.53 


Table 1. Solar element abund ance ado pted in this paper for calculating the 
chemical element abundance lAsnlund et al.l200td) . 


simulations have alrea dy been used in ICampisi et al] J 201 ll) and 
ISalvaterra et al] ]2013l) to study the properties of high -2 GRBs and 
their host galaxies. In the following, we will refer to the previous 
papers as C2011 and S2013, respectively. In the N-body hydro- 
dynamical chemistry calculations, the transition from a primordial 
(PopIII) to a PopII/I star formation regime is determined by the 
metallicity, Z , of the star forming gas. More specifically, a Salpeter 
initial mass function (IMF; oc m+ 2 ' 35 is the number frac¬ 

tion of stars per unit stellar mass) with mass range [100 — 500] Mg 
is adopted if Z < Z cr a = 10 _ 4 Zq, while a Salpeter IMF with 
mass range [0.1 — 100] Mg is used at larger metallicities. 

Besides gravity and hydrodynamics, the code follows early chem¬ 
istry evolution, cooling, star formation and metal spreading from 
all the phases of stellar evolution according to the suited metal- 
dependent stellar yields (for He, C, N, O, S, Si, Mg, Fe, etc.) 
and mass-dependent stellar lifetimes for both the pristine Po pIII 
regime and the metal-enriche d PopII/I regime ]Tomatore et al] 
l2007al : IMaio et al]|2007ll2010h . Yields for massive pair-instability 
supernovae (PISN) in the stellar mass range [140 — 260] Mg 
are taken by iHeger & Wooslevl J2002I ). Standard type-II super¬ 
nova (SNII) yields for progenitors with different metallicities 
and masses of [8 — 40] Mg are from IWooslev & Weaverl d 19951) . 
AGB metal production for lowe r-mas s sta rs is followed accord¬ 
ing to lvan den Hoek_&_Groenewegenl jl997i) and type-la SN (SNIa) 
yields are from Thielemann et al.ld2003l ). Several models for stellar 
nucleosynthesis are available in the literature and can give different 
yields; however, the overall metallic ity is usuall y affe cted in a mi¬ 
nor way (as already discussed in e.g. IMaio et alllolol) . 

In general, to better characterize the metal content of gas, stars or 
galaxies, it is convenient to define abundance ratios between two 
arbitrary species A and B as follows: 


[A/B]= 1 Ogl 0 (^)" 1 Ogl °(^)g 5 (1) 

where Na and Nb are the number densities of the two species and 
the sub script g refers to the corresponding solar content (Table [T| 
lAsplund et al.l2009l) . Abundance ratios are powerful tools to inves¬ 
tigate the pollution history and to probe stellar evolution models 
against observational data. 


3 RESULTS 

In this Section we will show results on the rate of GRBs originated 
from PopII/I stars exploding in a gas enriched by PopIII stars, to¬ 
gether with a discussion on the properties of the galaxies hosting 
such GRBs and observational strategies to identify them. This is 
an interesting issue, because high-redshift PopII/I GRBs exploding 
in a medium pre-enriched by former star formation episodes can 
reveal features of primordial PopIII generations through the metal¬ 
licity patterns of their hosting galaxy. 
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Figure 1. Redshift evolution of the ratio between the rate of GRB? and 
the total GRB rate. i=II—HII (black solid) and III (blue dotted-dashed for 
GRBIII up i and green dashed for GRBIII up 2 ). 


Figure 2. Rate of GRB? with redshift larger than 2 observable by 
Swift/BAT (see text for details). i=II (red dotted line), II—>111 (black 
solid) and m (blue dotted-dashed for (iR B111 1 and green dashed for 
GRBIIInj^)- 


3.1 Rate of GRBs 

We calculate the rate of GRBs as described in C2011. Here we 
provide only the essential information, while we refer the reader to 
the original paper for more details. The comoving GRB formation 
rate density, pGRBi, can be calculated as: 

pGRBi = fGRBi(,BH,ip* : i, (2) 


tgrbii^iii is as high as ~ 50% at 2 ~ 16, it becomes less than 
1% at 2 < 5. On the other hand, tgrbiii is always negligible, 
consistently with previous studies, as the transition from PopIII to 
PopII/I stars is very rapid. 

The (physical) observable rate of GRBi (in units of yr - 1 sr -1 ) at 
redshift larger than 2 can be written as (e.g. C 2011 ): 


where i indicates the subsample of GRBs considered, /GRBi is the 
fraction of BHs that can produce a GRB, C,bh,i is the fraction of 
BHs formed per unit stellar mass and p*,i is the comoving star for¬ 
mation rate density. We select star forming particles associated to 
PopII/I stars based on the criterium that Z > Z cr it , he. without an 
upper metallicity cut, as the GRB 1 in C2011. For the sake of clarity 
though, here we change the nomenclature and we refer to GRBII 
(GRBIII) as those GRBs originating from PopII/I (PopIII) stars. 
These are called GRBI (GRB3) in C2011. Following the original 
paper, we adopt C",bh,ii = 0.002 and / grbii = 0.028 for GRBII 
(which correspond to a minimum stellar mass of stars dying as BHs 
of 20 M@), and Cbh,iii = 0.0032 for GRBIII. fGRBin up2 = 
0.0116 is obtained assuming that no GRBIII is present among the 
266 GRBs with measured redshift, while the more stringent value 
fGRBin upl = 0.0039 is found imposing that none of the 773 
GRBs observed by Swift is a GRBIII. We additionally study a 
third class of GRBs, namely the GRBII that explode in a medium 
enriched by PopIII stars, which we indicate as GRBII—till. Here, 

C BH,II-*III = C,BH,II , fGRBI I — ill I = fGRBII , and p*,//->/// 
is explicitly calculated by selecting only those star forming parti¬ 
cles which have been enriched by PopIII stars. This task is accom¬ 
plished by using indicative metal ratios (see also next) for PopIII 
enrichment, consistently with the mentioned yield tables. In partic¬ 
ular, PopIII PISN enrichment results characterized by: [Fe/O] < 0, 
[Si/O] > 0, [S/O] > 0 and [C/O] < —0.5. These limits are used 
as reference selection criteria for PopII/I star forming regions en¬ 
riched by PopIII stars (see Sec. l3.2l for further discussion). 

Figure [T] shows the redshift evolution of the quantity tgrbi = 
pGRBi I (pGRBtot), where pGRBtot is the total GRB rate. As ex¬ 
pected, both tgrbiii and tgrbii^iii decrease with decreasing 
redshift, as the number of PopIII stars and their impact on the gas 
metal enrichment become less relevant. As a consequence, while 


Rgrbi{> z) = 


x 


lb 


L 


I. 


dz 




,dV(z') 1 PGRBi(z') 
dz’ 47 t (1 + z') 

dL'ipi(L'), 


(3) 


where 7 ?, = 5.5 x 1CP 3 is the beaming fraction for a n average jet 
opening angle of ~ 6 ° JGhirlanda et al J [20071 . l2013h . dV(z)/dz 
is the comoving volume element and ipi is the normalized GRB 
luminosity function. The factor (1 + z)^ 1 accounts for the cosmo¬ 
logical time dilation. The last integral gives the fraction of GRBs 
with isotropic equivalent peak luminosities above Lth, correspond¬ 
ing to an observed photon flux of ~ 0.4 ph s -1 cm -2 in the 15-150 
keV band of the Swift/BAT. We refer the reader to C2011 for a 
more exhau stive discussion on the choice of L t h and t/>; (see also 
Salvate rra e t SB. 

Interestingly, while there are much less GRBIII than GRBII—till 
(Fig.0, the former can be more easily detected because of their 
larger luminosity, as can be seen from Figure [2] resulting in 
very similar observable rates. Nevertheless, we should note that 
Rgrbiii has to be considered as an upper limit to the expected 
number of detections, as we have assumed here (as in C 2011 ) that 
all GRBIII can be detected by Swift/BAT thanks to their ex¬ 
pected extreme brightness. This should be true for GRBs origi¬ 
nating from massive PopIII stars (considered in this paper), while 
in case of more standard masses luminosities should be similar to 
those from regular PopII starsQ. 


1 GRBIII rates from C2011 have been updated by taking into account the 
new GRBs detected by Swift after the publication of C2011. These lead to 
a decrease of the C201 l’s PopIII GRB rate by a factor of 773/500 = 1.546 
for the ’upl’ case and by a factor of 266/140 = 1.9 for the 'up2’ case. 
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3.2 Metal signatures 

We now discuss the metal signatures which could help in identify¬ 
ing a GRBII—>111. Since the contribution to the metal enrichment 
from different stellar types is explicitly included in the simulations, 
the final abundance of each heavy element can be computed from 
the metal masses as traced at each output time. 

We use abundances normalized to oxygen, because this species 
is not very sensitive to theoretical yield uncertainties and is 
widely adopted in the literature. Additional reference species 
will be discussed in the following, though. For a given species 
X, we calculate the abundance ratio [X/O] in eq. |T| with 
Nx/No = (Aomx) / (Axmo), where Ax is the atomic mass 
number of element X and mx is the mass of X in a particle 
as given by the simulation. Moreover, it is possible to write 
log 10 (JVx/N 0 ) 0 = [log 10 (Nx/N h ) q - log 10 (No/N H ) e ]. 
Since GRBs are usually found in actively star forming environ¬ 
ments, only gas particles with non-zero SFR are considered while 
performing the calculations of the various chemical abundances. 
Before discussing the element abundances obtained from the sim¬ 
ulations, it is instructive to look at Figure[3] which shows the metal 
abundance of some elements from different stellar types of a simple 
stellar population, as expected by averaging over the corresponding 
IMF. More specifically, [X/O] in eq.[T]is obtained using Nx/No = 
[Ao / / [A x / Mo(m*)0(m*)dm*], 

where Mx is the stellar mass yield of element X for a given stellar 
population and the integration is performed over the relevant mass 
range for each stellar type. In the Figure, we do not show AGB 
stars, as their contribution to the metal enrichment is negligible 
for all shown elements, with the exception of C and O, which are 
produced in different amounts according to progenitor masses and 
metallicities, although typically [C/0]> 0. 

It is evident that each stellar phase is characterized by specific 
yields, that can be used to exclude or identify their contribution 
to metal enrichment. In particular, a positive and large value of 
[Fe/O] would indicate a strong contribution from SNIa, while 
a negative value could be due both to SNII and to PopIII stars 
(PISN). On the other hand, while SNII contribute with negative 
[Si/O] and [S/O], PopIII stars and SNIa contribute with positive 
ratios of the same elements. As mentioned above, AGB stars are 
the only ones to contribute with a usually positive [C/O] ratio. It is 
then clear that we can easily exclude gas enriched by SNIa as the 
one with [Fe/O]>0. Similarly, [S/O]<0 and [Si/O]<0 indicates 
gas enriched by SNII, and [C/O]>0 by AGB stars. It should be 
noted that the gas metallicity, Z, can affect the outflow of SNII 
and of AGB starfl but the considerations above still apply. So, to 
ensure that gas enriched by a stellar type other than PopIII stars 
is excluded, criteria such as [Fe/O]<0, [S/O]>0, [Si/O]>0 and 
[C/O]<0, could be adopted, although, as mentioned previously, 
the more stringent condition [C/O]<-0.5 has been chosen (see 
Sec.ITTI 

Figure [4] shows the distribution at various redshifts of some 
relevant elements in star forming particles with Z > Z cr it . With 
the exception of the highest redshift, when there are only a handful 
of star forming particles, two peaks can be clearly identified in 
each curve. The broader peak increases with decreasing redshift, 
and its metallicity abundance is similar to the one from SNII, so 

2 SNIa yields are weakly dependent on metallicities because they all 
derive from a white dwarf growing towards the Chandrasekhar mass 
and burning (in electro n-degenerate gas) C and O in equal proportions 
[Thiclc mann et al.l2004li . 



Figure 3. Element abundance as produced by: SNII (upper panel), PISN 
stars (middle) and SNIa (lower). The contribution from SNII depends on 
the gas metallicity, Z , as indicated in the legend. 


that metal enrichment of the particles in this peak might be led 
by gas expelled during SNII explosions. The increase of the peak 
with z is due to the time and metallicity evolution of the involved 
stellar systems. Indeed the relative shift in carbon abundance 
[C/O] is due to the metallicity-dependent carbon yields of massive 
stars augmented by the delay ed release of carbon from low- and 
intermediate-mass s tars (e .g. lAkerman et alj|2004l : ICescutti et alj 
120091 : ICooke et all 1201 ll) . Thus, the low values observed at 
2 ~ 17 — 13 are caused by SNII exploding in increasingly 
enriched gas (i.e. producing increasingly lower [C/O], as depicted 
in Fig. [3} and a subsequent rise at 2 ~ 13 — 5 led by AGB stars. 
The similar trend observed in [Fe/O] is a results of cosmic metal 
evolution as well, and at z ~ 5 — 6 the ratio is sustained also 
by the very first SNIa starting to explode after ~ 1 Gyr lifetime. 
Conversely, there is little variation in the distribution of the a ratios 
([Si/O], [S/O], [Mg/O]), as the nuclear reactions leading to their 
formation path are directly linked to O production. The element 
abundance of the narrow peaks seems more consistent with the 
yields from PopIII stars. Below z ~ 9 the peak stops increasing as 
a result of the extremely small PopIII SFR. 

Because some metallicity limits are equivalent in excluding 
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Element Abundance 

Figure 4. Abundance distribution (in number of particles) of, from top to 
bottom, C, Mg, Si, S and Fe of star forming particles with Z > Z cr a. The 
distributions are shown at z = 17 (solid red lines), 13 (green), 9.5 (blue), 
8.1 (black), 6.69 (magenta) and 5.25 (cyan). 


particles enriched by a specific stellar type (e.g. [Si/O], [S/O] and 
[Mg/O] can be equally used to exclude particles with metallicity 
dominated by SNII, while [C/O] or [Fe/O] can be used for SNIa), 
GRBII—HII can be identified also using less stringent criteria, 
and typically only two conditions are enough, such as [C/O] plus 
[Si/O], [C/O] plus [Mg/O], [Mg/O] plus [Fe/O] or [Si/O] plus 
[Fe/O], We find that if we adopt [C/O]<-0.5 and [Si/O]>0, or 



logloW©]) l°gi O (Z[Z 0 ]) 

Figure 5. Metallicity distribution of halos hosting a GRBII—tTII (solid 
black line) and of all PopII/I star forming halos (dotted red) at redshift 
z = 13 (top), 8.1 (middle) and 5.25 (bottom). Left panels: halos number 
fraction distribution normalized to the total number of star forming halos. 
Right panels: the distributions are “weighted” by their respective SFRs and 
normalized to the total SFR (see text for more details). The double arrows 
denote the Z range of GRB 050904 (solid black), GRB 130606A (dashed 
black) and GRB 111008A (dotted black). The black solid star denotes the 
Z of GRB 100219A. 


[C/O]<-0.5 and [Mg/O]>-0.4 the ratio rGRBii->m is ~ 1% 
(4%) and ~ 20% (50%) larger then the one shown in Figure [I] 
at 2 = 17 (5). At z > 12 the following criteria give similar 
results: [C/O]<-0.5 and [Mg/O]>-0.4, [Fe/O]<-0.1 and [Si/O]>0, 
[Fe/O]<-0.1 and [Mg/O]>-0.4. Finally, [S/O] and [Si/O] have 
similar selection effects, i.e. the same particles are selected if we 
use either [S/O]>0 or [Si/O] >0. 


3.3 Host galaxies 


In this Section we discuss some properties of galaxies hosting a 
GRBII—HII. Both observations and theoretical studies suggest that 
high-redshift (z > 5) long GRBs occur in objects with a_SFR lower 
than the one of the local Universe (z ~ Ol. lSalvaterra et alj (1201 30 . 
for example, find that such host galaxies have typically a stellar 
mass M* ~ 10 6 - 10 8 M 0 and a SFR ~ 0.003 - 0.3 M 0 yr' 1 . 
In general, metal enrichment from primordial star forming regions 
is efficient in polluting the structures hosting GRB episodes and in 
bringing typical metallicities above the critical value for stellar pop¬ 
ulation transition. This is in agreement with recent theoretical stud 


ies of primordial galax ies dMaio et alj|201 la 


imordial galaxie 

Wise et al Jl2012l 120141) and of GRB hosts dCampisi et alj 


a; Biffi_& Maio 

IlCampisi et al.1 


2013 


2011 


Salvaterra et al . 120131) . that suggest a typical metallicity range be¬ 
tween ~ 10 - 4 z 0 and ~ 10 1 Z 0 at z > 6. However, it is not 
clear yet how these statistical trends change when considering the 
sample of GRBII episodes in regions polluted by previous PoplII 
events. In order to discuss this issue, we compute the metallicity 
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distributions at different redshift for both the whole galaxy popula¬ 
tion and for the GRBII—>111 hosting population. 

In Figure [5] we show the metallicity distribution of halos hosting a 
GRBII—>111 (solid black line) and of all PopII/I star forming halos 
(dotted red line) at different redshift. While at the highest redshift 
both distributions are very similar in terms of halo number frac¬ 
tions (left panels), differences are visible at later times, when the 
typical metallicity of GRBII hosts increases. The simplest interpre¬ 
tation of these trends relies on the fact that primordial GRBII events 
take place in environments that have been previously enriched by 
the first PopIII generation and hence the two distributions almost 
coincide at early epochs, when only PopIII stars had time to ex¬ 
plode and pollute the gas. At later stages, PopIII stars play an in¬ 
creasingly minor role in the global star (and GRB) formation rate, 
which is dominated by PopII/I stars. This means that successive 
GRBII events will happen in gas that has been enriched by sev¬ 
eral generations of PopII/I stars, as already visible in the high -Z 
tail of the distributions at z — 8.1. Consequently, GRBII— >111 can 
survive only in the low-Z tail, where there is still some probabil¬ 
ity of having residual PopIII pollution in galaxies with weak star 
formation. This can be more clearly seen from the panels in the 
right columns, where we show a SFR “weighted” metallicity distri¬ 
bution, i.e. the fraction SFRi.j/SFR^tot, where SFR, -, is the SFR 
of halos in the j-th metallicity bin calculated for the component 
i =11, II—>111, and SFR^tot is the SFR of all halos for the same 
component. This explains the distributions at z = 5.25 and, more 
quantitatively, suggests a metallicity of roughly Z < lCT 2 ' 8 Z Q 
to highlight the presence of a GRBII—>111 (or possibly a GRBIII) 
in the early Universe. Another distinctive feature is the presence 
of two peaks at Z ~ 10 -4 Z 0 and IIP 2 ’ 8 Zq, found also in 
C2011 for GRBIII host galaxies (see their Fig. 3). This may be 
an indication that GRBII—>IIIs reside in galaxies similar to those 
hosting GRBIII, albeit with a higher star formation activity, as bet¬ 
ter highlighted by weighting the distribution of the host metal- 
licities by the SFR (right panels). Finally, we should note that 
a comparison w ith the metallicity range of the host galaxy of 
GRB 130606A l ICastro-Tirado et al.ll2013l: IChornock et al.ll20l 3l) 
seems to indicate that this is not a GRBII—>111. Although the red¬ 
shift of GRB 130606A is 5.91, i.e. slightly higher than what shown 
in the Figure, this does not change our conclusion because of the 
small redshift evolution of the peak at Z ~ 10 -2 — 10“ Zq. 
A similar conclusi on can be drawn for GRB 050904 at z — 6.3 
(Ka wai e t al . 20061) . GRB 111008A at z = 5.0 Isparre et al.ll2014l) 
and GRB 100219A at z = 4.7 dThone et al.ll2013l) . 

The corresponding SFR distributions, shown in Figure[6] span the 
range ~ IIP 4 ' 5 — 10 Mq yp 1 and neatly support the interpreta¬ 
tion, giving an upper limit of ~ 10~ 2 Mq yr _1 for GRBII—>111 
host galaxies (higher values of the SFR are found, but with a much 
lower probability). More specifically, the SFR displays a trend sim¬ 
ilar to the one of the metallicity, with GRBII—>111 forming prefer¬ 
entially in halos with a SFR smaller than the one of GRBII hosts, 
and a large fraction of halos with SFR ~ 10 -4 — 10 -3 Mq yr _1 
hosting a GRBII—>111 at all redshifts. The highest probability of 
detecting a GRBII—>111, though, is in halos with a slightly higher 
SFR, i.e. ~ 10 -2 ' 5 Mq yr , as shown in the right column (note 
that here the distributions are calculated as in Fig. [5] where Z bins 
are substituted by SFR bins). 

Finally, in Figure [7] we show the stellar mass distribution of ha¬ 
los hosting a GRBII—>111 (solid black line) and of all PopII/I star 
forming halos (dotted red) at different redshift. Despite both dis¬ 
tributions being very similar in terms of halo number fraction at 
the highest redshift (left panels) as z decreases GRBII—>III are still 



log 10 (SFR [M 0 yr- 1 ]) log 10 (SFR [M 0 yr" 1 ]) 

Figure 6. SFR distribution of halos hosting a GRBII— >-111 (solid black line) 
and of all PopII/I star forming halos (dotted red) at redshift z = 13 (top), 
8.1 (middle) and 5.25 (bottom). Left panels: halos number fraction distribu¬ 
tion normalized to the total number of star forming halos. Right panels: the 
distributions are “weighted” by their respective SFRs and normalized to the 
total SFR (see text for more detail). 



log 10 (M, [Mq]) log lo (M„[M 0 ]) 


Figure 7. Stellar mass distribution of halos hosting a GRBII— >-111 (solid 
black line) and of all PopII/I star forming halos (dotted red line) at redshift 
z = 13 (top), 8.1 (middle) and 5.25 (bottom). Left panels: halos number 
fraction distribution normalized to the total number of star forming halos. 
Right panels: the distributions are “weighted” by their respective SFRs and 
normalized to the total SFR (see text for more details). 
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preferentially hosted by halos with M* ~ 10 5 while the typ¬ 
ical host of a GRB has a stellar mass about one or two orders of 
magnitude larger. The right panels of the Figure, showing the dis¬ 
tributions of the host stellar masses “weighted” by their respective 
halo SFR (note that here the distributions are calculated as in Fig. [5] 
where Z bins are substituted by M* bins), confirm that the high¬ 
est probability of detecting a GRBII—► III is, as above, in structures 
with M* ~ 10 5 M 0 at all redshifts. 


4 CONCLUSIONS 


We use N-body hydrodynamical chemistry cosmological simula¬ 
tions to compute the rate of high-redshift GRBs arising from the 
death of PopII stars and exploding in an environment previously 
enriched by PopIII supernovae (GRBII—>111). Indeed, the detection 
of these objects can provide an indirect but powerful probe of the 
existence and the properties of the first stars in the Universe. 

The results presented in this work are based on numerical cal¬ 
culations including detailed chemistry evolution, gas cooling ac¬ 
cording to suited metal-dependent stellar yields (for He, C, N, O, 
S, Si, Mg, Fe, etc.) and mass-dependent stellar lifetimes for both 
the pristine PopIII regime and the metal-enrich ed PopII/I regime 
dTornatore et alj l2007al : iMaio et all 120071 . l201Ch . Metal pollution 
proceeds from dense star forming regions to surrounding lower- 
density environments and is responsible for the enrichment of the 
Universe already in the first billion years, when early massive stars 
die and possibly generate GRB e vents . Numer ical limitations (see 
a more complete discussion in e.g. lMaio et al J201cL and references 
therein) in addressing such issues might lie in the uncertainties of 
metal-diffusion modeling (both in SPH and grid-base approaches). 
Effects due to the PopII/I IMF slope and features or to the particu¬ 
lar critical metallicity adopted in the range ~ 1CP 6 — 10 -3 Z 0 are 
minor. 

On the other hand, a PopIII IMF shifted t owards low er masses 
as recently suggested by some authors (Yoshida etjdJ_ 20071 : 


ISuda & Fuiimotol2010l ; lci ark et alj20lU ; IStacv & Bromml2014l ) is 

expected to have a larger impact on our results. In this case, in fact, 
the metal yields from PopIII stars are expected to be more similar 
to those from PopII/I stars and thus it would be increasingly diffi¬ 
cult (if at all possible) to discriminate gas enriched by primordial 
stars based on metal abundances. We explicitly checked this issue 
by re-running the same box with a PopIII IMF in the 0.1 — 1OOM 0 
mass range and we found results consistent with analogous studies 
based on DLA obse rvations and employing low-mass pristine IMFs 
(e.g. iKulkami et aill2014l) . More specifically, we found that while 
assuming a SN mass range of 8-40Mq it would be extremely dif¬ 
ficult to distinguish GRBII—HII because the metal yields of PopIII 
stars are too similar to those of PopII/I stars, including stars as 
massive as 1OQM 0 would induce ex tremely high Oxygen and Car¬ 
bon yields l lHeger & Wooslevll2010l) . In the latter case, GRBII—>III 
could be identified e.g. using the condition [Si/O]<-0.6, [S/O]<-0.6 
and [C/O]>-0.4. 

Possibly different metal yields as predicted by different stellar 
structure models could slightly shift the resulting abundance ratios 
according to the specific input yield tables, but no m ajor changes 
are expected, mostly at low metallicities dCescutti et al.l2009l) . Res¬ 
olution effects are generally unavoidable, but they are likely to 
play a minor role for the abundance ratios explored in this work 
and for the star formation regime transition, since these processes 
are predominantly determined by stellar evolution timescales, 
which are constrained relatively well. The cosmology adopted is 


ACDM, but changes in the background cosmological scenario and 
the introduction of high-order non-linear effects would not af- 


feet dramatically the baryonic properties within bound structures 

(see e.g. Maio et alj 2006|, 201 lbl. 20li Maio & Iannuzzi 

201 ll: 

Maio & Khochfar 2012;; de Souza et al. 2013; Maio & Viel 

2014). 


Despite that, some particular dark-matter models might cause a de- 
lay in early structure formation (as is the case for warm dark matter; 
Ide Souza et al]|2013l; IMaio & Viell2014l) . 

We have found that the rate of GRBII—*111 rapidly decreases with 
cosmic time from ~ 50% at 2 = 16 down to less than 1% at 2 < 5. 
We expect that one tenth of the GRBs exploding at 2 = 10 should 
report the metal signature of the enrichment by a PopIII supernova 
explosion. By convolving the intrinsic rate of these events with 
the GRB luminosity function, we found that ~ 0.06 GRBII—*111 
yr -1 sr _1 should be sufficiently bright to trigger Swift /BAT. 
This means that ~ 0.8 GRBII—till should be present in the en¬ 
tire Swift database. 

The identification of such events can be obtained by looking at pe¬ 
culiar metal abundance ratios. We showed that [C/O] and [Si/O] 
alone could be enough to distinguish GRBII—HII from other GRB 
populations. However, in practice, the detection of more elements 
(such as S and Fe) could help observers in further confirming the 
nature of the source. While we have presented all our results with 
respect to abundance ratios normalized to oxygen, because O is not 
very sensitive to theoretical yield uncertainties, it is possible to get 
particularly interesting hints by referring metal abundances to Si, as 
well. In this case, PopII star forming particles enriched by PopIII 
stars can be identified as those with [C/Si] < —0.5, [O/Si] < 0 
and [Fe/Si] < —0.4, where the first two conditions exclude gas 
enriched by SNII, while the last one excludes pollution from SNIa. 
AGB stars produce no Si and thus they are excluded in any case 
from such analysis. We find that the above conditions give results 
very similar to those obtained using our reference criteria in the en¬ 
tire redshift range. 

It should be noticed that it is hard to get accurate Fe abundances 
from observatio ns at high redshift, while C, Si and O are usu ally 
easier to probe JCastro-Tirado et al]l2013l ; [Ka wai et ai] l2006l) . In 
addition , Fe yiel ds m ight be overestimated of up to a factor of ~ 2 
dCescutti et ailj 120091 ). Thus, Fe is not the best suited element to 
identify GRBII-HIIs. 

We also explored the properties of galaxies hosting GRBII—HII 
in our simulations. These are very similar to th ose that have a 
high probability to host a GBRIII dSalvaterra et alj 120131) . i.e. 
GRBII—HII are typically found in galaxies with stellar mass (~ 
10 4 ' 5 - 1O 7 M 0 ), SFR (~ 1(T 3 - O.lM 0 yr _1 ) and metallicity 
(~ 10 -4 — 1O -2 Z 0 ) lower than those of galaxies hosting a GR¬ 
BII. It is worth to note that galaxies with log 10 (Z/Zq) < —2.8 
most probably host GRBII—HII and GRBIII. A GRB with mea¬ 
sured metallicity below this Z limit should be considered as a 
strong candidate for being originated in a PopIII enriched environ¬ 
ment (or by a PopIII progenitor). 

From an observational point of view, GRBII—HII signatures can 
be directly investigated for the currently known high-2 GRBs. At 
2 > 6 metal absorption lines in the optical-NIR spectrum have 
been detected and metal abundance ratios meas ured only for two 
bursts, namely GRB 05090 4 at 2 = 6.3 dKawai et al .1120061) and 
GRB 130606A at 2 = 5.91 dCastro-Tirado et alj|2013h . 

The observed abundances for GRB 050904 are [C/H]=-2.4, 
[0/H]=-2.3, [Si/H]=-2.6 and [S/H]=-1.0 translating into [C/0]=- 
0.1, [Si/O]=-0.3 and [S/O]=0.7. On this basis, we can ex clude 
GRB 050904 to be a GRBII— > III. Recent studies by [Thone et alj 
d2013h have reported an updated abundance for S of [S/H]=—1.6± 
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0.3 and a corresponding metallicity of log(Z/Zp) = — 1.6 ±0.1, 
instead of log(Z/Z©) =: —1.3±0.3 ilCawai et al . 20061) as inferred 
from Si absorption lines of the afterglow spectrum. However, such 
variations do not af fect our conclusion on GRB 050904. 

For GRB 130606A ICastro-Tirado et all i2013l) reported a 3a up¬ 
per limit of [S/H] < —0.82 and determined the lower limits 
of the oxygen and silicon abundance as [Si/H] > —1.80 and 
[O/H] > —2.06, and metall icity Z as ~ (1/60 — 1/7) Zq. 
Results from Ichomock et a (2013) give [Si/H] > —1.7 and 
[S/H] < —0.5. while lHartoog et al l2014h find [C/H] > —1.29, 
[O/H] > -1.88, [Si/H] = -1.33 ± 0.08, [S/H] < -0.63 and 
[Fe/H] = —2.12 ± 0.08. Using these information, it is possible 
to estimate corresponding limits of [S/O] < 1.24, [Si/O] < 0.55, 
and [Fe/O] < —0.24. These are too weak, though, to probe the 
nature of the environment in which the burst explodes. We also 
note that th e lower lim its on the total metallicity based on Si and O 
dChornock et a 1 ][2~0 1 3l) suggest a moderately enriched environment 
where GRBII are likely to dominate. However, the possibility that 
GRB 130606A is associated to a PopIII enriched environment can 
not be completely excluded, although a comparison with the metal¬ 
licity expected for a galaxy hosting a GRBII—HII seems to indicate 
otherwise. 

Anoth er possible candidat e at z = 4.7 could be GRB 100219A, for 
which [Thone et all ( t2013h reported the following measurements: 
[C/H] = —2.0±0.2, [O/H] = —0.9±0.5, [Si/H] = -1.4T0.3, 
[S/H] = —1.1 ± 0.2 and [Fe/H] = —1.9 ± 0.2. Metallicity is 
about 0.1 Zq and is derived from the abundance of sulphur. From 
such values one can get [C/O] = —1.1±0.7, [Si/O] = —0.5±0.8, 
[S/O] = —0.2 ± 0.7 and [Fe/O] = —1 ± 0.7. The precision at 
la-level is not enough to draw robust conclusions, though. 

We also note that preliminary analyses are available for 
GRB 111008A at z — 5.0 JSparre et all l2014ll . which is the 
highest-redshift GRB with precise measurements of metallicity, 
column densities and iron, nickel, silicon and carbon fine-structure 
emissions. Reported abundance ratios within 2a-error bars are: 
[Si/H] > -1.97, [S/H] = —1.70±0.10, [Cr/H] = -1.76±0.11, 
[Mn/H] = -2.01 ± 0.10, [Fe/H] = -1.74 ± 0.08, [Ni/H] = 
— 1.64 ± 0.19, [Zn/H] = —1.58 ± 0.21. The authors do not give 
carbon and oxygen values, but they note that [Fe/H], [Ni/H] and 
[Si/H] measurements suggest metallicities around ~ 1 — 6 per cent 
solar, well above the 10 2 8 Zq limit for GRBIII and GRBII—HII. 
Morevover, the inferred reddening (Av = 0.11 ±0.04 mag ) of the 
afterglow indicates that the dust-to-metal ratio of the host galaxy is 
0.57±0.26, similar to typical values of the Local Group. Hence, the 
GRB environment should be strongly dominated by PopII/I enrich¬ 
ment, while it is unlikely that it has been pre-enriched by massive 
PopIII stars. 

Finally, the metallicity inferred from X-ra y measures at higher z 
dCampana et al .1 [20 111 : Starling e t alj|2013h suggests values above 
a few percent solar. These might not be reliable though, as the X- 
ray value of the H column density they are based on is likely af¬ 
fected by int ervening metals along the line-of-sight dCampana et alj 
I2012L l2015ll. Wit h future X-ray facilities like the Athena satel- 
lite jNandra et ahl l2013l Fl it will become possible to directly mea¬ 
sure the abundance patterns for a variety of ions (e.g. S, Si, Fe) 
superimposed on bright X-ray afterglows. This, in principle, will 
allow to discriminate between different nucleo-synthesis sources 
djonker et all2013h and identify GRBII-HIIs. 


3 http://www.the-athena-x-ray-observatory.eu/ 
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